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The microscopic structure of neutron-rich W-shell nuclei is investigated by using the spherical shell-model 
in the sd-pf valence space with the extended pairing plus quadrupole-quadrupole forces accompanied by the 
monopole interaction (EPQQM). The calculation reproduces systematically the known energy levels for even- 
even and odd-mass nuclei including the recent data for ''■'S, and Ar. In particular, the erosion of the N = 2S 
shell closure in "^^Si can be explained. Our EPQQM results are compared with other shell model calculations 
with the SDPF-NR and SDPF-U effective interactions. 
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I. INTRODUCTION 

The current experimental and theoretical investigation fo- 
cuses on the study of evolution of nuclear structure around 
the shell gap at N ~ 28. The N = 2% shell closure is a tradi- 
tional one in the nuclear single-particle spectrum driven by the 
spin-orbit interaction. For example, the '^^Ca nucleus has been 
known to exhibit double-magicity as a result of the neutron 
subshell gap separating the /7/2-orbit and the /5/2P-shells at 
N = 2%. However, recent theoretical studies and experimental 
data have questioned the persistence of the traditional magic 
numbers, and revealed that the = 28 shell gap is eroded 
when moving away from the stability line. The experimental 
data for the "^^S and '^^Si isotones indicate a clear breaking of 
the N = 2% magicity because of the observed small 2+ en- 
ergy 1 1-3]. Recently, a direct evidence of collapse of the = 
28 shell closure due to the level inversion between the lowest 
7/2^ and 3/2^ levels has been observed in "^^^S |4]. In addition, 
toward — 28, the degeneracy between the lowest 1/2+ and 
3/2+ states in the odd-proton P, CI, and K isotopes suggests 
15] that the shell gap between the 1^3/2 and 51/2 proton orbits 
almost collapses at = 28. The reduction of the N = 2% neu- 
tron gap combined with this degeneracy is regarded as the pos- 
sible origin of erosion of the N = 2% magicity in the neutron- 
rich iii-shell nuclei. 

To gain the insight for the evolution of shell structure in 
the neutron-rich icZ-shell nuclei and to understand the details 
of collapse of the N = 2E magicity, shell-model calculations 
in the full sd-pf space are desirable. However, conventional 
shell-model calculations in the full sd-pf space are not pos- 
sible at present because of the huge dimension in the con- 
figuration space. To study the neutron-rich .sii-shell nuclei, 
therefore, one needs to incorporate truncations in the valence 
space. One choice is such that valence protons are restricted 
to the sd shell and neutrons to the sd-pf shell with twelve sd 
frozen neutrons. Such an approach was oijiginally introduced 
with the effective interaction SDPF-NR The nuclei 

in this mass region are known to have a variety of shapes 
(e.g. shape coexistence discussed in Refs. |2, 4]), and the nu- 
cleus ^^Si has been the subject of debate on its magic nature 
lUKTll]. The SDPF-NR interaction was applied to the neutron- 



rich S(i-shell nuclei, and one found that it describes well the 
excitation energies for very neutron-rich isotopes with Z > 14. 
However, the calculation deviated from the experimental data 
for the silicon (Z = 14) isotopes where the calculated ener- 
gies for ^^^^^Si were too high. To improve the results, the 
pZ-shell effective interaction in SDPF-NR was renormalized 
to cornpensate for the absence of 2p-2h excitations from the 
core II12I1 . With the reduction of pairing in the renormaliza- 
tion for Z < 14, the 2^ excitation energies in silicon isotopes 
agreed nicely with experiment. Thus, two SDPF-U interac- 
tions, one for Z > 14 and the other for Z < 14, are respec- 
tively needed for the description of collectivity at A^ = 28 in 
the isotopic chains of sulfur (Z =16) and silicon (Z =14). In a 
similar way, the SDPF-NR2 and SDPF-NR3 interactions have 
recently been suggested to describe the sudden change in nu- 
clear structure at A^ = 28 ifTsll . From the above discussion, one 
sees that a consistent shell-model treatment for collectivity at 
the N = 2S shell closure in different isotopic chains has not 
been successful so far. 

The main purpose of the present article is to construct a 
unified effective interaction for the neutron-rich ^li-shell nu- 
clei, which can be consistently applicable to both Z > 14 
and Z < 14 isotopic chains, and to understand the shell evo- 
lution in this exotic mass region. It is well known that re- 
alistic effective interactions used in the low-energy nuclear 
structure study are dominated by the pairing plus quadrupole- 
quadrupole (P + QQ) forces with inclusion of the monopole 
term lfl4ll . As documented in the literature, the extended 
P+QQ model combined with the monopole interaction works 
well for a wide range of nuclei 1 15, 16]. This effective inter- 
action is called hereafter EPQQM to distinguish it from shell 
models with other effective interactions. The EPQQM model 
has demonstrated its capability of describing the microscopic 
structure in different N ~Z nuclei, as for instance, in the fp- 
shell region |15] and the /pg-shell region |16]. Recently, it 
has been shown that the EPQQM model is also applicable to 
the neutron-rich Cr isotopes II 1711 . 

The monopole interaction plays an important role in our 
discussion. The monopole shifts in the spherical shell model 
have been introduced to account for the non-conventional 
shell evolution in neutron-rich nuclei. Connection between 
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the monopole shifts and the tensor force flS*] has been studied 
within the self-consistent mean-field model using the Gogny 
force ifigll . It has recently been clarified that the possible phys- 
ical origin of these interactions is attributed to the central and 
tensor forces Il20l I2TI1 . For example, when one goes from 
the Ca down to the Si isotopes a significant reduction of the 
N = 2Si shell gap can be found. By removing protons from the 
t/3/2 orbit in nuclei between Ca and Si, the strong attractive 
proton-neutron monopole force between the ndi, and v/7 /2 
orbits is no longer present 1^^,122]. This makes the neutron 
fi 12 less bound, and thus reduces the size of the = 28 gap 
in nuclei between '^''Ca and ^'^Si. 

The paper is arranged as follows. In SecUll we outline our 
model. In Section |III1 we perform the numerical calculations 
and discuss the results for the neutron-rich nuclei in the sd-p f 
shell region. Finally, conclusions are drawn in Section lTV] 
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FIG. 2: (Color online) Experimental and calculated first excited 2^ 
energies E{2^) for the neutron-rich A' = 28 isotones. 



II. THE MODEL 

We start with the following form of Hamiltonian, which 
consists of pairing and quadrupole-quadrupole terms with the 
monopole interaction 



H — //sp + //j 



' Hp 



Hqq + i/m 



J=Q.2 ' 



Mk 



M 

+ EL^m(fl^) L^/Mr/f(«^M/Mr/f(fl^), (1) 

a<h T JMK 

where b in the third term is the length parameter of harmonic 
oscillator. We take the J — Q and 7 = 2 forces in the pairing 
channel, and the quadrupole-quadrupole iQQ) forces in the 
particle-hole channel [15, 16 ] . The monopole interaction is 
denoted by //m^ where the global monopole force is neglected 
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FIG. 1: (Color online) Comparison of experimental and calculated 
energy levels for the odd-mass nuclei near the neutron or proton shell 
closure in the sd region. 




because it does not affect the excitation energies of the low- 
lying states. This isospin-invariant Hamiltonian ([TJ is diago- 
nalized in a chosen model space based on the spherical basis. 
We employ the shell-model code ANTOINE |23] for the nu- 
merical calculation. In the present work, we consider an '^O 
core and employ the sd-pf model space comprising the Qd^ 12, 
1^1/2' 0'^3/2 active proton orbitals and the 0d^i2, l«i/2. ^dj,/!, 
l/'3/2' 0/5/2' active neutron orbitals with twelve 
sd frozen neutrons. We employ the same single-particle en- 
ergies as those of Ref. |12]: 6^/5/2 = 0.0, £51/2 — 0.78, 
e,/3/2 = 5.60, e/7/2 = 9.92, £,,3/2 = 10.01, £^1/2 = 15.15, and 
S/5/2 — 10.18 (all in MeV). We adopt the following interac- 
tion strengths for the EPQQM forces: go = 0.540,^2 = 0.678, 
and Xi — 0.474 (all in MeV). Finally, for the monopole terms, 
the strengths are chosen to be (all in MeV): 

kl='{fi/2,d5/2) = ~03, 

klr^{d5/2,P3/2) = -0.2, 

'lr° {^1/2, P3/2) = -0.6, 

^m=H/7/2,/5/2) = -0.1, 



k^= 


"°(/7/2 7^^3/2) 


= -0.65, 


k^= 


'° ('^5/2,/' 1/2) 


= 0.55, 


k^= 


"°('^3/2j/'3/2) 


= 0.5, k 


k^= 


'^{ds/2,dil2) 


= -0.2, 


k^= 


'\f7/2,Pl/2) 


= -0.01, 


k^= 


"'(/7/2J./7/2) 


= -0.3, 



{fy2,Py2) = -0.03, 
{P3/2,P3/2) = -0.45.(2) 



The rest of the monopole terms are neglected in the present 
calculations. The above EPQQM force strengths are deter- 
mined so as to reproduce the experimental energy levels for 
odd-mass nuclei that have one particle or one hole on top 
of the neutron and proton shell closures. The results are 
compared with experimental data |24] in Fig. [T] Here the 
root mean square (RMS) deviations between the experimen- 
tal and theoretical excitation energies are 0.24, 0.14, and 0.20 
(in MeV) for the 3/2^, 1/2^, and 5/2^ states, respectively, 
which means that the agreement is good. We emphasize that 
in order to describe the data correctly, a very important part of 
the interaction is those monopole forces. The attractive T — 
monopole interactions 7td^/2'^fi/2 and 7^d5/2'^f7/2 have re- 
cently been proposed to be attributed to the monopole effect 
of the central and tensor forces in the nucleon-nucleon inter- 
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FIG. 3: (Color online) Effective neutron single-particle energies at 
A' = 28 from Z = 8 to Z = 20. 
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FIG. 4: (Color online) Energy splitting between the lowest 7/2 and 
3/2^ states in neutron-rich N = 21 isotones I4l l25ll27h . 



action 1120112 lH . where the central term produces a global con- 
tribution and the tensor term generates local variations. They 
are responsible for the collapse of the N — 22i shell closure in 
the neutron-rich iii-shell nuclei, which gives rise to the drastic 
change in nuclear shape of the N = 2E isotones. In fact, it has 
been reported that the tensor force produces the pronounced 
oblate minimum in the potential energy surface for "'^^Si [221 . 
Along the K isotopic chains, the repulsive T^d^/i-Vp^/i and 
attractive ?r^i/2-V/:>3/2 monopole interactions lead to an inver- 
sion of the expected ordering in the lowest 3/2+ and 1/2+ 
levels in ^^K. In addition, these interactions reproduce the evo- 
lution of the lowest 3 /2^ state in ^^Si. The monopole interac- 
tions between nd^ and the /7/-shells act so as to reproduce 
the low-lying states in "^^Ca. 



III. NUMERICAL RESULTS AND DISCUSSIONS 

A. The shell erosion in A' = 28 

We discuss changes in the shell structure when moving 
away from the valley of stability. Large energy gaps at the tra- 



ditional magic numbers in stable nuclei may be washed out in 
the neutron- or proton-rich regions. The fundamental question 
of how the nuclear shell structure evolves with proton or neu- 
tron excess is one of the main motivations to study nuclei far 
from stability. The study of shapes in such exotic nuclei serves 
as a sensitive test for the predictive power of nuclear models. 
Figure |2] shows the systematics of the first excited 2+ states 
for the N — 2% isotones, which are compared with three shell 
model calculations, respectively with the EPQQM, SDPF-U, 
and SDPF-NR effective interactions. As one can see, all the 
calculations are capable of producing the drastic drop in en- 
ergy when removing protons from the double magic nucleus 
48Ca. 

To understand the erosion of the N = 2% shell gap, the neu- 
tron effective single-particle energies (ESPE) as a function of 
proton number are shown in Fig. [3] The v/7/2 and V/73/2 
orbitals at Z = 20 becomes degenerate at Z = 16. This de- 
generacy is due to the attractive proton-neutron monopole in- 
teraction between the v/7/2 and nd^i2 orbitals. As seen in 
Fig. m the inversion of the lowest 7/2^ and 3/2^ states in ^^S 
(Z = 1 6 and = 27) just reflects the degeneracy of the two or- 
bitals, v/7 /2 and v/73/2, at Z = 16 in Fig. [3] On the other hand, 
the neutron ESPE in the SDPF-U interaction does not show 
such a strong reduction of the N = 2% shell gap from Z = 20 to 
Z = 16. At Z = 16 the gap stifl equals about 3.5 MeV, and the 
inversion between the 7/2^ and 3 /2^ states is ascribed to the 
combined effect of the gap reduction due to the monopole in- 
teractions and the increase of the multipole correlations. The 
N = 2% shell closure is eroded in ''^''Ar and '^^S, after the re- 
moval of only two and four protons, respectively. This rapid 
disappearance of rigidity in the N — 2% isotones has been as- 
cribed to a reduction of the neutron shell gap at = 28 com- 
bined with that of the proton subshell gap at Z = 16, leading to 
increased probability of quadrupole excitations within the fp 
and sd shells for neutrons and protons, respectively. In fact, 
the occurrence of the low-lying isomer at 320.5 keV in ^^S 
has been interpreted in the shell-model framework as result- 
ing from the inversion between (v7/2^)^' and (v3/2^)+' 
configurations ll25i l26ll . We have thus understood that the in- 
version has the origin of vanishing N = 2% magicity at Z = 16. 

In Fig. |5] the proton ESPE with Z 20 are shown as a func- 
tion of neutron number It is seen that the shell gap between 
the Kdi,i2 and 7r5i/2 proton orbitals in the double closed-shell 
nucleus ^''Ca (Z — N = 20) decreases with increasing neutron 
number. The 7zdj,/2 orbit is almost degenerate with the 7r.si/2 
orbit atN = 28, and lies just below the nsi /2 orbit at A^ = 28. 
An inversion of the 7zd^/2 and 7r5i/2 orbits thus occurs above 
A^ = 32. The effect of adding neutrons to the 77/2 orbital is 
primarily to reduce the proton ^ 1/2-1^3/2 gap. The relevant or- 
bitals, 7tsi/2 and Kd^i2, are known to become degenerate at 
A^ = 28. This degeneracy enhances the quadrupole correlation 
energy of the configurations with open neutron orbitals. The 
ESPE from A^ = 34 to A^ = 40 for the nd^ /2 are different from 
those of the SDPF-U calculations. It would be attributed to 
the monopole effect between the 7td^i2 and the vfsi2- Fig- 
ure|6]shows the energy splitting between the lowest 3 /2+ and 
1 /2+ states as a function of mass number for the neutron-rich 
isotopes of phosphorus, chlorine, and potassium. The calcu- 
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lations are in a good agreement with available experimental 
data [23, for the phosphorus and chlorine isotopes, but 
not for the potassium isotopes. The energy splitting between 
3/2+ and 1/2+ for K isotopes in Fig. |6] deviates from the 
experimental data. This would be related with the disagree- 
ment of the first excited 2+ enegies in Fig. |2] The variation 
trend reflects the reduction of the ndi,i2-7isii2 splitting in Fig. 
|5]as neutron number increases. The inversion of the 11(1^12- 
%sii2 orbits at 28 corresponds to the observation that the 
ground-state of ^^K is actually the 1/2+ state, not the expected 
3/2+ from the hole state (%d^i2Y^ ■ The near degeneracy of 
Tid^j2 and ns\j2 is again attributed to the attractive proton- 
neutron monopole interaction between the v/7/2 and %dT,i2- 
The erosion of the N = 2% shell gap is enhanced at Z = 16 by 
the degeneracy of the 7^^3/2 and 7r.Si/2 proton orbits. In Fig. 
|6] we also show the energy splitting between the lowest 3/2+ 
and 1 /2+ states for the "-45^ jm^lei L28J. 
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FIG. 5: (Color online) Effective proton single-particle energies at 
Z = 20 from = 20 to A' = 40. 

C. Level scheme and B{E2) values 



B. Comparison with other shell model calculations 



A consistent description of collectivity at = 28 for dif- 
ferent isotopic chains of the neutron-rich sd-shell nuclei has 
been a challenge for shell model calculations. The early cal- 
culations with the SDPF-NR interaction described the level 
scheme and transitions for these isotopic chains. However, 
it was found that the calculated results deviate from the ex- 
perimental data. For the silicon isotopes the calculated en- 
ergy levels are too high. To improve the agreement, Nowacki 
and Poves have proposed two effective interactions [12], one 
(SDPF-U) for Z < 14 and the other (SDPF-Ul) for Z > 14. 

It is interesting to compare the energy levels obtained from 
the shell model calculations with different effective interac- 
tions (EPQQM, SDPF-U, SDPF-Ul, and SDPF-NR). Figure 
|7] shows the first excited 2+ energies as a function of neu- 
tron number for the Ca, Ar, S, and Si isotopes. The EPQQM 
calculations agree well with the known experimental data, ex- 
cept for ^'^Si and ^'^S. The SDPF-NR calculations also agree 
with the data for the Ca and Ar isotopes, but cannot reproduce 
those of the Si isotopes. The calculated first excited 2+ states 
lie higher than experiment. On the other hand, the SDPF-U 
(SDPF-Ul) interaction fails to reproduce the known experi- 
mental data for Z > 14 (Z < 14). Moreover, considerable dif- 
ferences between the EPQQM and SDPF-NR (and also SDPF- 
U, SDPF-Ul) can be seen around the neutron number = 34, 
where no experimental data are currently available to discrim- 
inate the predictions. We note that none of the interactions 
can reproduce the first excited 2+ energy level of ^'^Si with 

= 20. This suggests that it would be necessary to include the 
neutron excitations from the icZ-shell to the p f -shell across the 
N = 20 energy gap. The disappearance of the N = 20 magic 
structure between the sd and pf shells has already been estab- 
lished for the neutron-rich nuclei with N = 20 isotones. 



In this subsection, we present the theoretical level schemes 
and electromagnetic E2 transition probabilities calculated 
with the EPQQM interaction, and compare the results with 
experiment. In Fig. [8] we show energy levels for the even- 
even sulfur isotopes, where for each isotope, the left and right 
part of energy levels denote the experimental data lfl3il24l[30ll 
and the calculated energies, respectively. The agreement be- 
tween the calculation and experiment is fairly good. In partic- 
ular, the calculated first excited 2+ energy level for ^^S repro- 
duces very well the recently observed data [13]. Very recently, 
energy spectrum for ^^Ar has been experimentally obtained 
lUlS. In "^^S, the calculated £'(4+)/£'(2+) ratio is about 
2.5, which qualitatively suggests a transitional or 7-soft nature 
for this nucleus. In '^^S, the £'(4j^)/£'(2j^) ratio is 3.0, which 
is quite close to the rigid rotor value. The situation signifi- 
cantly differs in '^'^S where both deformed and spherical con- 
figurations are predicted to coexist and mix weakly with each 
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FIG. 6: (Color online) Energy splitting between the lowest 3/2+ and 
1/2+ states in neutron-rich P, CI, and K isotopes. 
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FIG. 7: (Color online) Comparison between experimental and calculated first excited 2+ energies £{2'^) with several effective interactions, 
EPQQM, SDPF-U, SDPF-Ul, and SDPF-NR, for Si (Z = 14), S (Z = 16), Ar (Z = 18), and Ca (Z = 20) isotopes. 
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FIG. 8: (Color online)Experimental and calculated energy levels for 
even-even Sulfur isotpes [2., 24 , 29 , 30.1. 



Other, which is supported by the observation of the low-lying 
Oj^ state at 1.326 MeV f30]. In the shell model calculation, 
this state lies below the 2^ state. The presence of a low-lying 
state is considered as a signature of a spherical-deformed 
shape-coexistence and the present data support the weakening 
of the N — 22i shell gap. Figure |9] shows the energy levels for 
the even-even argon isotopes. One finds that the calculated 
energy levels can reproduce the experimental data ll24l [33ll . 
The first excited 2+ energy E{2^) for Ar isotopes decreases 
with increasing neutron number from ^^Ar to ^'^Ar, but with 
a sudden increase for the N — 2S isotope. This suggests the 
persistence of the = 28 gap in Ar nuclei. In Figs. |8]and 
ill we note that the RMS deviations for the 2^, 2j, and 4| 
states are respectively 0.22, 0.71, and 0.92 (in MeV), and the 



agreement becomes worse with increasing excitation energy. 
The disagreement may require the further improvement of the 
EPQQM interaction. 

To demonstrate the EPQQM model calculation for odd- 
mass nuclei in this mass region, we show in Fig. [TO] energy 
levels for some odd-mass S and Ar isotopes. In the compari- 
son, the experimental energy levels for "^^S and '^^Ar are taken 
from the recent experiments P.'sT]. From the calculation, the 
ground state of ^^Sis3/2^, not the expected 7/2 ^ , suggesting 
that this nucleus is deformed. The occurrence of a low-lying 
isomer at 320.5 keV in ^^S is interpreted in the shell-model 
framework as resulting from the inversion between the natu- 
ral (v/7/2) ' and the intruder configurations. Both 
the low B(£'2;7/2^ — > 3/2^) value and the absence of cal- 
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FIG. 9: (Color online)Experimental and calculated energy levels for 
even-even Argon isotopes L24, . .3 1.-33.1 . 
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TABLE I: Calculated S(£2 : 2+ 
Ar, and Si isotopes. 



0+) values for the even-even S, 



FIG. 10: (Color online)Comparison of calculated e nergy le vels with 
the recent experimental data for S and Ar nuclei fll M H ^IsljSTll . 



culated deformed structure built on the 7/2^ isomer suggest 
a coexistence of different shapes in the low-lying structure of 
the nucleus ^^S. In ^^Ar, the ground state 3/2^ and the ex- 
cited state 1/2^ can be interpreted as a coupling of a pt,/2 and 
Pi/2 neutron to the 0+ ground state of ^^Ar. The 5/2^ and 
l/2~ levels result from the coupling of a neutron to the 
2+ state of 46 Ar. 

In Table I, we compare B{E2) values between the theoret- 
ical calculations and experiments for the S, Ar, and Si iso- 
topes. The experimental B{E2) values are taken from Refs. 
iH m HI HHi. In our B{E2) calculations, we take the ef- 
fective charges e^i for protons and ey for neutrons as follows; 
— 1.35e and gy = 0.35e for S isotopes, e„ = 1.50e and 
= 0.50e for Ar isotopes, and e„ ~ l.l5e and ey = 0.15e for 
Si isotopes. On the other hand, in the SDPF-NR and SDPF-U 
calculations the effective charges e^t — 1.35e and gy = 0.35e 
suggested in Ref. II2I1 are used for S and Si isotopes, and 
the effective charges e^r = 1 -SOe and ey = 0.50e are adopted 
for Ar isotopes. The B{E2) values were previously calculated 
and discussed with the SDPF-NR interaction by Retamosa et 
al. 101 • As one can see, our EPQQM calculation is in a good 
agreement with the experimental B{E2) values for the ^^^"^^S, 
but larger than the experimental value for ^^s. On the other 
hand, both the SDPF-U and SDPF-NR calculations also re- 
produce the data B{E2) except for ^^ 445 jjj Xable I, we also 
compare the calculated B{E2) values with data for the Ar iso- 
topes. The experimental B{E2) values are taken from Refs. 
HmmiM. One observes that the experimental B(£'2) in- 
creases monotonically from Ar to the midshell nucleus ^^Ar, 
and then decreases towards ^^Ar. We note that all the theoret- 
ical calculations agree fairly well for ^^-^'^Kx^i for 46 Ar the 
calculated B{E2) values are all too large ll32ll4l|] . All the cur- 
rent shell model calculations with different interactions yield 
a similar value, which may indicate that this is a robust result. 
Very recently, however, a new experimental value 1 14 eHrn^ 
has been reported f^^, and is in a good agreement with all 
the shell-model predictions. In Table I, this new experimental 
value is included for comparison. It is apparent from the re- 
duction in excitation energy of the 2'\ states and the enhance- 
ment in B{E2) that collective features are gradually developed 
in neutron-rich isotopes between N = 20 and N = 28. How- 
ever, a different approach, the self-consistent mean-field cal- 
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48Ca 




5 


3 


3 



culation by Werner et al. [43], yields a much smaller B{E2), 
in a disagreement with the new experimental data. We com- 
pare B{E2) values between theoretical calculations and exper- 
iment for the Si isotopes. The agreements with the available 
data are satisfactory. The B{E2) values increase with increas- 
ing neutron number A^. The calculated B{E2) values are quite 
similar for the EPQQM and SDPF-U interactions. For SDPF- 
NR, however, the predicted values are smaller for the heavier 
isotopes 40-42si^ 

Finally in Table II, we present calculated spectroscopic 
quadrupole moments for the N = 2% isotones. For all the ef- 
fective interactions, small, positive 2's are found for 4**Ca, 
indicating a spherical character of this nucleus. For the 44s 
nucleus, its small 2'^ energy, large B{E2) value, and the pres- 
ence of a second excited isomer at low excitation suggest a 
mixed ground state configuration with spherical and deformed 
shapes. Recent experiments on the odd mass sulfur isotopes 
around 44s seem to indicate that 44s is a deformed nucleus 
with strong shape coexistence whereas 42s can be considered 
as a well-deformed system. The enhancement in B{E2) of 
46 Ar coincides with a shape change from a small prolate de- 
formation in 44 Ar to a large deformation in 46 Ar. The pre- 
dicted spectroscopic quadrupole moment Qs{2'l) in 44 Ar is 
somewhat larger than the experimental value, but with a cor- 
rect sign. In 42si, the yrast sequence of Oj'", 2'^, does not 
follow the rotational 7(7+ 1) law. However, its quadrupole 



7 



moments are consistent with a deformed oblate structure. The 
ground state of '^''Mg and that of ^^Si are predicted to be 
prolately and oblately deformed, respectively. All the shell 
model calculations indicate consistently a rapid variation in 
quadrupole moment for the = 28 isotones, suggesting a 
drastic shape change in this mass region. 

IV. CONCLUSIONS 

We have investigated the microscopic structure of the 
neutron-rich sd-shell nuclei using the spherical shell model 
with a new effective interaction in the sd-pf valence space. 
The interaction is of the extended pairing plus quadrupole- 
quadrupole type forces with inclusion of the monopole inter- 
action. The calculation has reproduced reasonably well the 
known energy levels for even-even and odd-mass nuclei in- 
cluding the recent data for ^^S, ^^S and '^^Ar A special atten- 
tion has been paid to the N = 2S shell-gap erosion. It has been 



shown that the attractive T = monopole interactions Ttdj, /2' 
v/7/2 and 7^df,i2-vfi 12 very important for the variation of 
the shell gap at = 28. This shell evolution is considered 
as the monopole effect of the central and tensor forces in the 
nucleon-nucleon interaction. The strong attractive monopole 
interaction l^dT^j2-V ft 12 causes the inversion of the 3/2+ and 
1/2+ states in odd-proton nucleus '^^K. However, the inver- 
sion is recovered in the neighboring nucleus '^^K, which is due 
to the repulsive monopole interaction 7idi,i2-Vp^i2- We have 
studied the shape change in the neutron-rich ==28 isotones, 
and predicted the drastic shape change when protons are re- 
moved away from the Z = 20 isotone ^^Ca, i.e. going towards 
the neutron drip-line. 
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